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Abstract 



We study the stochastic resonance phenomenon in the overdamped two cou- 
pled anharmonic oscillators with Gaussian noise and driven by different exter- 
nal periodic forces. We consider (i) sine, (ii) square, (iii) symmetric saw-tooth, 
(iv) asymmetric saw-tooth, (v) modulus of sine and (vi) rectified sinusoidal 
forces. The external periodic forces and Gaussian noise term are added to one 
of the two state variables of the system. The effect of each force is studied 
separately. In the absence of noise term, when the amplitude / of the applied 
periodic force is varied cross- well motion is realized above a critical value (/c) 
of /. This is found for all the forces except the modulus of sine and recti- 
fied sinusoidal forces. For fixed values of angular frequency u of the periodic 
forces, /c is minimum for square wave and maximum for asymmetric saw- 
tooth wave, /c is found to scale as ^e'^'^^'^ + B where A and B are constants. 
Stochastic resonance is observed in the presence of noise and periodic forces. 
The effect of different forces is compared. The stochastic resonance behaviour 
is quantized using power spectrum, signal-to-noise ratio, mean residence time 
and distribution of normalized residence times. The logarithmic plot of mean 
residence time tmr against 1/{D — Dc) where D is the intensity of the noise 
and Dc is the value of D at which cross-well motion is initiated shows a sharp 
knee-like structure for all the forces. Signal-to-noise ratio is found to be max- 
imum at the noise intensity D = D-^^ at which mean residence time is half 
of the period of the driving force for the forces such as sine, square, symmet- 
ric saw-tooth and asymmetric saw-tooth waves. With modulus of sine wave 
and rectified sine wave, the SNR peaks at a value of D for which sum of 
Tmr in two wells of the potential of the system is half of the period of the 
driving force. For the chosen values of / and a;, signal-to-noise ratio is found 
to be maximum for square wave while it is minimum for modulus of sine and 
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rectified sinusoidal waves. The values of at which cross-well behaviour is 
initiated and -Dmax are found to depend on the shape of the periodic forces. 

PACS: 02.50.-r; 05.40.-a; 05.45.-a 

Keywords: Over damped two coupled anharmonic oscillators; stochastic res- 
onance; periodic forces; power spectrum, mean residence time; signal-to-noise 
ratio. 
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Typeset using REVT^ 



1. Introduction 

The phenomenon of stochastic resonance has simulated series of theoretical, numerical 
and experimental works on the dynamics of multi-stable systems subjected to both periodic 
and noisy driving. In many nonlinear systems, the occurrence of stochastic resonance has 
been studied with the external periodic force being of the form f sin cut [1-16]. There are few 
studies with aperiodic forces. For example, enhancement of response of certain nonlinear 
systems by the added noise with aperiodic forces is observed [17-20]. Jamming (undesirable 
noise) occurring in an array of FitzHugh-Nagumo oscillators was found to be suppressed via 
stochastic resonance by the input aperiodic signals such as weak amplitude modulated, fre- 
quency modulated and chaotic signals [21]. Stochastic resonance can happen in the absence 
of periodic excitation [5, 22] and noise [23, 24] also. The former was called spontaneous 
stochastic resonance and latter as deterministic stochastic resonance. 

In recent years, the dynamics of nonlinear systems with different periodic forces has been 
studied [25-30]. Particularly, the onset of homoclinic chaos in a damped pendulum driven by 
a periodic string of pulses modulated by Jacobian elliptic function and periodic 5-function is 
investigated [25] . The existence of a new multiple-period-doubling bifurcation route to chaos 
and modified bifurcation structures in the Duffing oscillator perturbed by periodic pulses are 
reported [26]. Anti-control of chaos is accomplished by adding certain periodic forces [27]. 
A nonlinear feedback controller that generates chaotic behaviour in an oscillator driven by 
a distorted force has been proposed [28]. Suppression of chaos by 5-pulse in Duffing- van der 
Pol oscillator [29] and by a variable shape pulse function in a coupled pendulum-harmonic 
oscillator system [30] is found. 

Certain dynamics are studied with forces other than / sin cut. It is important to study a 
particular dynamics with different types of forces and make a comparative study of effects 
induced by them. This is because different forms of periodic forces can be generated and 
easily applied to real mechanical systems and electronic circuits. Further, the study of 
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influence of various forces in nonlinear systems will be helpful to choose a suitable force in 
creating and controlling nonlinear behaviours. 

Motivated by the above, in the present paper, we wish to study stochastic resonance 
phenomenon with different periodic forces in the overdamped two coupled anharmonic os- 
cillators given by [31] 

X ^ aix -bix^ + Sxy'^ + F{t) +rj{t), (la) 
y = a2y- b2y^ + 5x^y. (lb) 

In eqs(l) F{t) and ri{t) represent the external forcing term and the noise term respectively. 
In the absence of external periodic force, noise and damping terms the potential of the two 
coupled anharmonic oscillators is 

y (X, y) = - + '^y^ - ^xV- (2) 

In a very recent work, Baxter et al [32] considered the eqs(l) as a model for competition 
between two species. They analysed the probability of survival of the species using the path 
integral method in eqs(l) without F{t) and adding noise to both the state varaibles x and y. 
Different types of intermittent-lag synchronization are shown to occur in a parametrically 
driven two coupled anharmonic oscillators [33]. 

The organization of the paper is as follows. In section 2, to be self-content we give the 
mathematical form of the different external forces used in our study. The effect of periodic 
forces without noise term is discussed in section 3. For fixed values of the parameters, the 
system is found to exhibit cross-well periodic orbit above a certain critical value (/c) of the 
amplitude of the forces. The value of /c depends on the shape of the forces, fc is found to 
increase with increase in the frequency cu of the forces and scales as Ae^^'^^'^ + B. In section 
4 we consider the system in the presence of external periodic force and noise. For each force 
except the modulus of sine and rectified sinusoidal wave stochastic resonance is observed 
at Dmax at which mean residence time is half of the period of the driving force. Stochastic 
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resonance is characterized using power spectrum, signal-to-noise ratio, mean residence time 
and distribution of normalized residence time. We compare the effect of various forces on 
stochastic resonance phenomenon in terms of its characteristic quantities. Section 5 contains 
summary and conclusion of our present study. 

2. Types of periodic forces 

The periodic forces of our interest are 

1. sine wave 

2. square wave 

3. symmetric saw-tooth 

4. asymmetric saw-tooth 

5. modulus of sine wave 

6. rectified sine wave 

Figure (1) shows the shape of the above forces. Throughout our study, we fix the period Tq 
of all the forces as 2t:/u!. 

The mathematical representations of the periodic forces are the following: 

(1) Sine wave 

The sine wave is given by the form 

F^^{t)^f sin cut, (3) 

where / is the amplitude of the force and cu is the angular frequency of the force. 

(2) Square wave 

The mathematical form of square wave is 



/, (2n - 2)7i/uj <t<{2n- 1)ti/uj 

(4) 

-/, {2n- 1)'k/u <t <2n'K/uj, n = l,2,... 



Fsq{t) = 

For numerical implementation, we rewrite the force (4) as 
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Fsqit) = 



/, <t < n/u 



(5) 



-/, 7l/u) < t < 2tt/ijj, 
where t is taken as mod(27r/c<;). That is, Fsq(t) = Fsq{t + 2tt/u). Its Fourier series is 

(3) Symmetric saw-tooth wave 

Symmetric saw-tooth wave is mathematically represented as 



(6) 



^sstW=<l 



f, {2n-2)n/u;<t<{^)n/u; 
-f + 2/, {^)n/u;<t<{^)n/u; 



(7) 



4n- 



To ' V 2 



-) n/u <t < 2n7r/u;, n = 1, 2, . . . 



where Tq is the period of the force. The numerical implementation of the force is given by 



^sst(^) = 



^, 0<t< 7r/2, 



w 



^ + 2/, 7r/2u; < i < 37r/2u; 



(8) 



^ - 4/, 37r/2a; < t < 2t^/u, 
where t is taken as mod(27r/a;). Its Fourier series is 

(_l)(n+i) sin(2n- iVi 



^sstW = 4E 



n=l 



(2n - 1)2 



(9) 



(4) Asymmetric saw-tooth wave 

The mathematical representation of the wave form shown in fig. Id is 



^ast(^) 



{2n -2)t: / uj <t <{2n-l)T:/ uj 
^-2f, {2n-l)n/u <t <2n7r/LJ, n = l,2, 



(10) 



The force (10) is rewritten as 

^ast(^) 



if' 0<i<7r/ 



To 

Ml 

To 



UJ 



¥ - 2/, 7r/u;<t< 27r/a;, 



(11) 



where t is taken as mod(27r/a;). The Fourier series of i^ast(^) given by 



^ n=l ^ 

(5) Modulus of sine wave 

The modulus of sine wave is given as 

^msiW = /|sinM/2)|. (13) 

Its Fourier series is 

2/ 4/^ cosn^t 
Fms.it) ---- }^^J^^;^y (14) 



(6) Rectified sine wave 

The mathematical representation of rectified sine wave (fig. If) is 

/, (2n-2)n/uj <t < (2n- l)n/u} 

Frsiit)-{ (15) 
I 0, {2n-l)n/uj <t <2mT/u;, n = l,2,... 

The numerical implementation of the force (15) is 

/, <t < tt/cu 



^rsi(^) = <! 

0, n/u <t < 2^/00. 
The Fourier series of rectified sine wave is 



(16) 



^ , , 1 2 J^^ cos 2//L^'f siiiLL-'t , 
for < t < n/u and for n/u <t< 2n/u. 

3. Cross- well motion due to different periodic forces without noise term 

In this section we considered the system (1) in the absence of noise term rj{t) but with 
different external forces. The parameters of the system are fixed at ai = 1.0, 02 = 1.1, 
bi — 1.0, 62 — 1.0. For this choice, the potential is a four-well potential. The potential wells 

are designated as for x > 0, y > 0; y+_ for x > 0, y < 0; V-+ for x < 0, y > 0; V for 

X < 0, y < 0. The period of the all the forces is set to 27r/a;. For each force the response of 
the system is studied by varying the amplitude of it. 



First, we consider the effect of the sinusoidal force fsinut. For a fixed value of oj and for 
small values of /, four period- Tq {27i/u) orbits one in each potential well exist. The size of 
the orbit increases with increase in the value of /. At a critical value of /, /c, the trajectory 
makes a cross-well motion forming a coupled orbit. However, the period of the cross-well 
orbit remains as Tq. Since the external forcing is added only to x-component of the system, 

the coupled orbit traverse the wells V"++ and V-+ or the wells V+- and V depending upon 

the initial conditions. For / in the range [0, 50], chaotic behaviour is not observed. This 
is confirmed from the bifurcation diagram and Lyapunov exponent calculations. Similar 
behaviour is observed for various values of cu in the interval [0, 1] with step size cu = 0.05 
and also for other forces except for modulus of sine wave and rectified sine wave. In the 

case of modulus of sine wave added to x-component of the system the wells y_+ and V 

alone are oscillating. The wells V++ and V+_ are stationary. For initial conditions in the 

well V+- or V , after the transient evolution the system remains in the well V+-. In the 

case of rectified sine wave during the first half of the period of the drive cycle only the wells 

V_+ and V oscillate. Therefore, cross-well behaviour is not observed for / in the range [0, 

50] in both the cases. 

For the force f sin cut with lu — 0.05, the critical value of / (at which cross- well motion 
is initiated) is numerically found to be 0.409. For forces such as square, symmetric saw- 
tooth and asymmetric saw-tooth waves /c is 0.393, 0.4408 and 0.475 respectively. Figure 2 
shows the variation of fc with lu. In this figure the dots represent numerical result and the 
continuous lines are the best curve fit. We find that fc increases with increase in cu. For a 
fixed uj the fc for square wave is minimum while it is maximum for asymmetric saw-tooth 
wave. That is, cross-well motion is initiated at a lower value of / for square wave when 
compared with other forces. The cross-well motion occurs relatively at a higher value of / 
in the case of sine, symmetric saw-tooth and asymmetric saw-tooth waves. For all the forces 
fc is found to vary with cu as Ae^''^^'^ + B. The values of (^4, B) obtained for square, sine, 
symmetric saw-tooth and asymmetric saw-tooth forces are (0.471, —0.130), (0.754, —0.409), 
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(0.987, -0.602) and (1.284, -0.834) respectively. 

4. Stochastic resonance with different periodic forces 

In the previous section, we studied the effect of different periodic forces separately in 
the absence of noise term r){t). Now, we include the noise term and study the stochastic 
resonance phenomenon by varying the noise intensity D. Here again we perform the analysis 
for each force separately. The noise term is chosen as Gaussian random numbers. The 
parameters are fixed at Oi = 1.0, 02 = 1.1, 61 = 1.0, 62 = 1-0, S — 0.01 and u; — 0.05. In 
the absence of noise, the critical value of / at which cross-well orbit occurs for square, sine, 
symmetric saw-tooth and asymmetric saw-tooth waves are 0.393, 0.409, 0.4408 and 0.475 
respectively. The value of / is fixed below fc- We fix / = 0.38 for all the forces so that in 
the absence of noise the motion is confined to single well alone, that is, there is no cross- 
well motion. We integrated eqs(l) using a fourth-order Runge-Kutta method with step size 
At — {2t^/uj)/N, N — 2001. Noise is added to the state variable x as Xj+i — > Xi^i-\-\fKiDC,{t) 
after each integration step /St. Here C(t) represents Gaussian random numbers with zero 
mean and unit variance and D denotes noise intensity. 

4.1. Cross-well motion and mean residence time 

First, we illustrate the effect of noise in the presence of square wave. Cross- well motion is 
not realized for D < 0.0002. Ai D — D^, orbit switching from one-well to another well 

is initiated. For example, cross-well behaviour is observed between the two wells and 
a.t D = Dc- Similar behaviour is observed when the square wave is replaced by other 
wave forms with same cu and / values. However, the value of Dc is found to be different 
for other forces. D^. values for square, sine, symmetric saw-tooth, asymmetric saw-tooth, 
modulus of sine and rectified sine waves are 0.0002, 0.0005, 0.002, 0.0025, 0.05 and 0.03 
respectively. D^. is relatively low for the square wave. 

For each periodic force, the switching between the wells become very rare for D values 
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just above Dc- However, as D increases switching also increases. To characterize this 
behaviour we numerically calculated the mean residence time. Residence time tr is defined 
as the duration of the trajectory of the system residing in a well (for example, before 
jumping to another well (say, and vice-versa. Mean residence time tmr in the wells 

y_+ and are computed for a set of 10^ residence times. Rapid variation of tmr is found 
when the noise intensity D is varied from D^. Figure 3a shows In(TMR) versus 1/{D — Dc) 
when the applied force is square wave. Later, we will show that signal-to-noise ratio (SNR) 
is maximum at this value of D. Different form of variation of tmr is found for D < i^max 
and D > L'max- The curve in the fig. 3a has a sharp knee ai D = L>max- At this value 
of D the derivative of tmr is discontinuous. Sharp-knee shape like variation of maximal 
Lyapunov exponent versus control parameter is observed for band- merging crisis [34] . Mean 
residence time tmr for the other periodic forces are also calculated. The logarithmic plot of 
Tmr versus 1/{D — Dc) is shown in figs.3b-d for sine, symmetric saw-tooth and asymmetric 
saw-tooth forces. All the plots resemble the knee-like structure. The values of -Dmax are 
found to be 0.14, 0.1, 0.12 and 0.17 for square, sine, symmetric saw-tooth and asymmetric 
saw-tooth forces. The maximum SNR values are obtained at these values of D for the above 
forces. 

For the forces depicted in figs.la-d the sign of the forces is positive during one half of 
the cycle and negative during the second half of the cycle. For the forces in fig.le and If the 
sign of the force is never negative. This difference in the forces and the symmetry between 
the wells have infiuence on the mean residence time. For the forces shown in figs.la-d for 
D values just above -Dmax the tmr in the wells and are same and greater than the 
period Tq — 2-k/uj ~ 125.67 of the drive cycle. As D increases tmr decreases. In the case of 
modulus of sine and rectified sine waves, the mean residence time tmr in the well K-+ for 
D just above Dc is less than To/2 = tt/o; « 62.83 while it is greater than Tq for the other 
well That is, tmr is not same for the two wells V^_+ and V++. tmr in both the wells 

decreases with further increase in D. Figures 4 and 5 show the variation of tmr with D for 
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the modulus of sine and rectified sine waves respectively. 

4.2. Stochastic resonance 

Time series plot in the presence of external square wave for few values of noise intensity 
D is shown in fig.6. For small values of D the system exhibit the same periodic behaviour of 
noise free case but slightly perturbed by the noise. For D < D^ — 0.0002, the trajectory of 
the system resides in a single well as in fig.6a. At the critical value Dc the trajectory jumps 
randomly from one well to another. In fig. 6b for D — 0.001 just above the state variable 
X switches irregularly and rarely between the wells and In the presence of forcing, 
the system initially in the well, say, V++ is forced by the noise to leave the well. Then the 
system enters the well and wanders irregularly there for some time and jumps back to 
the well V++ and so on. For D = 0.001 tmr is computed as fa 790. The switching is not 
periodic. The residence times are randomly distributed, tmr decreases with increase in D. 
For D — 0.14, Tmr is found to be ^ tt/o; ^ 62.83. In this case nearly periodic switching 
between the wells y++ and y_+ occur. There is a co-operation between the periodic driving 
force and the noise. The x-component of the trajectory switches between the positive and 
negative values with the period approximately half of the period of the applied external 
periodic force. This is clearly seen in fig.6c. 

For large values of D, the motion is dominated by the noise term. Figure 6d shows 
the trajectory for D = 2.5 where the trajectory jumps erratically between the wells. That 
is, loss of coherence is produced by the large noise intensity. For D in the range [D^, 4], 
jumping motion between the wells and V-+ alone observed. The system (1) is studied 
with different periodic forces as a function of noise intensity. We observed similar response 
of the system for all the forces except modulus of sine and rectified sine waves. 

The trajectory of the system exhibited periodic transitions between the wells at L>max = 
0.1, 0.12, 0.17 for sine, symmetric saw-tooth and asymmetric saw-tooth waves. Synchro- 
nization of the input wave form with the output at these values of D is shown in fig. 7 for 
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various forces. In this figure for all the forces ai D — £)max we can notice a common feature. 
Nearly at the end of one half of a drive cycle the trajectory in one well is likely to jump to 
the another well and after the next half cycle it is likely to return back. This is the signature 
of stochastic resonance. 

4.3. Chctracterization of stochastic resonance by signal-to-noise ratio 

We characterize the stochastic resonance phenomenon in terms of SNR and the peak 
value at the half driving frequency in the distribution of normalized residence times. SNR 
is calculated from the power spectrum of the state variable x. To compute power spectrum, 
we have used fast Fourier transform technique [35, 36]. The power spectrum is obtained 
using a set of 2^° data collected at a time interval of (27r/a;)/10. To get more accurate power 
spectrum we computed it for 25 different realization of Gaussian random numbers and then 
obtained the average spectrum. 

In the absence of noise term, the power spectrum of x component of the system with 
different periodic forces is first analyzed. With square wave as the external forcing, the 
power spectrum has peaks at few odd and even integral multiples of the forcing frequency. 
When the noise is included, the peaks at odd integral multiples are alone dominant. With 
increase in noise intensity D, the height of the peaks increase up to D — Dmax (at which 
tmr ~ tt/o;) and then decrese. The power spectrum for four values of D in the presence of 
square wave is shown in fig. 8. Dominant peaks are seen at fl = 0.05 and 0.15. 

Signal-to-noise ratio is calculated from power spectrum using the formula 

SNR = 101ogio(^/iV) dB. (18) 

In eq(18) S and N are the amplitudes of the signal peak and the noise background respec- 
tively. S is read directly from the power spectrum at the frequency cu of the driving periodic 
force. To calculate the background of the power spectrum about cu we considered the power 
spectrum in the interval Q, — uu — Auj and uj + Auu after subtracting the point namely Q, — cu 
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representing the stochastic resonance spike. The average value of the power spectrum in 
the above interval is taken as background noise level at to. We have chosen Ao; = 0.01. 
Figure 9a shows the plot of SNR as a function of D for square wave. The calculated SNR 
increases rapidly with the noise intensity D, peaks at the critical value Dmax = 0.14 and 
then decreases slowly for D > -Dmax- With modulus of sine and rectified sine waves for 
small values of noise intensity D, rare switching between the wells and V-+ is observed. 
However, periodic switching with the half of the period of the drive cycle is not observed 
when D is varied. Figure 9b shows SNR versus D for various applied forces. With the 
forces such as sine, square, symmetric saw-tooth and asymmetric saw-tooth waves, SNR 
has a peak at D — Dmax at which tmr ~ tt/o; and there is almost a periodic switching 
between the wells y_+ and V^^ as shown in fig. 7. For the forces given in figs.le-f the tmr in 
is always less than tt/lu. However, stochastic resonance is observed at the value of D at 
which the sum of the mean residence times of both the wells V-+ and V++ is n/u;. This 
happens ai D — 0.35 and D — 0.3 for modulus of sine and rectified sine waves respectively. 
For these forces SNR initially decreases with increase in D from Df. for a while, increases 
with increase in D, reaches a maximum value and then decreases with further increase in 
D. SNR is found to be highest for square wave. Modulus of sine wave and rectified sine 
wave yields the minimum SNR value. 

4.4. Characterization of stochastic resonance by probability distribution of nor- 
malized residence times 

The probability distribution of normalized residence times denoted as P can also be used 
to quantify stochastic resonance. Normalized residence time is calculated as follows. For a 
fixed noise intensity D the residence time tr in a well is computed for a set of 10^ transitions. 
Then, normalized residence times are obtained by dividing tr by the period Tq of the applied 
force. Then we calculated P{tr/Tq). P is found to be same for the wells V_+ and VL^. In 
fig. 10 we have plotted P versus tr/Tq for four values of noise intensity with the applied 
periodic force being square wave. We observe a series of Gaussian like peaks, centered at 
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odd integral multiples of To/2. The heights of these peaks decrease exponentially with their 
order n. The strength Pi of the peak at tr = To/2 is a measure of the synchronization 
between the periodic forcing and the switching between the wells. The variation of Pi can 
be used to characterize the stochastic resonance phenomenon. For values of D nearly above 
Dc, tr/Tq is distributed relatively over wide interval of time (fig. 10a). As D increases the 
range of tr decreases and P of smaller tr increases. At P* = Pmax as shown in fig. 10c the 
noise induced oscillations have period To/2 and hence P{tr/Tq = 1/2) becomes maximum. 
As D is further increased Pi decreases. Similar behaviour is observed for the periodic forces 
given in figs.la-d. Figure 11 shows Pi versus D for various forces. The effect of various 
forces on Pi can be clearly seen in this figure. For small values of D, tr/Tq is distributed 
over wide interval of time. Therefore, Pi increases with increase in D, reaches a maximum 
Sit D = Z^max and then decreases. With the modulus of sine and rectified sine waves, the 
peak Pi at tr/Tq = 1/2 is numerically calculated for various values of D in the wells V-+ 
and For the well the value of Pi decreases with increase in the value of D. But 
for the well Pi increases with increase in D reaches a maximum at which SNR is 

maximum and then decreases. Figures 12 and 13 show the plot of Pi versus noise intensity 
D for the forces in fig.le and fig. If respectively. 

5. Summciry and conclusion 

We have studied numerically the effect of different periodic forces and Gaussian noise in 
the overdamped two coupled anharmonic oscillators eqs(l). In order to compare the effect 
of external forces we have fixed the period of the forces as 27r/a;. The average of all the 
forces over one drive cycle is zero except for modulus of sine wave and rectified sine wave. 
The system shows periodic behaviour in the presence of external periodic forces and in the 
absence of noise term. Since the forcing is added only to x-component, the system performs 
transitions between the wells y++ and y_+ only. The value of /c is found to be small for 
square wave when compared with other forces. Cross-well behaviour is delayed in the case 
of sine, symmetric saw-tooth and asymmetric saw-tooth forces. With the other two forces, 
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cross-well motion is not observed. Stochastic resonance is observed with all the forces used. 
The mean residence time, power spectrum, SNR and distribution of normalized residence 
times are used to characterize stochastic resonance. The features of stochastic resonance are 
found to be same for the forces such as sine, square, symmetric saw-tooth and asymmetric 
saw-tooth waves. Maximum SNR is obtained at the value of D at which sum of the residence 
times ~ tt/o; for modulus of sine and rectified sine waves respectively. The values of and 
-Dmax are found to be different for different periodic forces. Signal-to-noise ratio values are 
found to depend on the type of forcing used. Square wave gives maximum SNR. Modulus 
of sine wave and rectified sine wave yields minimum SNR. The height of the peak Pi at 
-Dmax obtained from normalized residence time distribution vary with the periodic forces. 
In the nonlinear dynamics literature, many interesting results have been obtained with the 
quasiperiodic, chaotic, aperiodic forces and with the different noise terms. Therefore, it is 
interesting to study the stochastic resonance behaviour with the above forces. 
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FIG. 1. Form of different periodic forces: (a) sine wave, (b) square wave, (c) symmetric 
saw-tooth wave, (d) asymmetric saw-tooth wave, (e) modulus of sine wave and (f) rectified sine 
wave. For all the forces the period Tq is 27r/a; with u = 0.05 and the amplitude is set to 1. 
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for various forces. The curves (a) 
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sine wave, symmetric saw-tooth wave and asymmetric saw-tooth wave respectively. Dots represent 
numerical data and continuous lines are the best curve fits. 
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FIG. 3. Mean residence time in logarithmic scale versus 1/{D — Dc) for (a) square wave, (b) 
sine wave, (c) symmetric saw-tooth wave and (d) asymmetric saw-tooth wave. The parameters of 
the system are ai = 1.0, 02 = 1.1, h = 1.0, 62 = 1.0, 6 = 0.01, / = 0.38 and uj = 0.05. The symbol 
X indicates the value of l/(-Dmax — -Dc) where -Dmax is the value of D at which tmr = tt/o;. 
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FIG. 6. Time series plot of eqs(l) in the presence of square wave for few values of noise intensity 
D. The parameters of the system are ai = 1.0, 02 = 1.1, 61 = 1.0, 62 = 1-0, S = 0.01, / = 0.38 and 
u = 0.05. The applied square wave is also plotted. 
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FIG. 7. Synchronization of the output signal with the different applied periodic forces such as 
(a) sine wave, (b) square wave, (c) symmetric saw-tooth wave and (d) asymmetric saw-tooth wave. 
In all the subplots the value of D is Dmax- 
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FIG. 8. Power spectral density plot in the presence of square wave for few values of noise 
intensity D. 




FIG. 9. (a) Signal-to-noise ratio plot for range of values of noise intensity D in the presence 
of square wave, (b) Signal-to-noise ratio plot for various periodic forces such as sine (continuous 
curve), square (dotted curve), symmetric saw-tooth (— •— curve), asymmetric saw-tooth (dashed 
curve), modulus of sine (dots joined by a continuous curve) and rectified sine (circles joined by a 
continuous curve) used. In all the cases / = 0.38 and a; = 0.05. 
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FIG. 10. Normalized residence time distribution for few values of D in the presence of square 
wave with amplitude / = 0.38 and frequency a; = 0.05. 
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FIG. 11. Height of the peak Pi in the normalized residence time distribution at tr = To/2 for 
square wave (marked by •), sine wave (o), symmetric saw-tooth wave (*) and asymmetric saw-tooth 
wave (x). 
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FIG. 12. Pi versus D for the wells (a) F_+ and (b) with the applied force as modulus of 
sine wave. 
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FIG. 13. Pi versus D for the wells (a) F_+ and (b) with the applied force as rectified sine 



wave. 
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